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ABSTRACT 

The  Pareto  distribution  has  been  proposed  recently  as  a  model  for  intensity 
clutter  measurements,  for  maritime  high  resolution  radar  returns.  Using  the 
theory  of  spherically  invariant  random  processes,  the  Neyman-Pearson  optimal 
detector  is  derived.  The  generalised  likelihood  ratio  test  yields  a  simple  sub- 
optimal  approximation,  and  its  performance  is  gauged  against  the  whitening 
matched  filter. 
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Coherent  Multilook  Radar  Detection  for  Targets  in  Pareto 

Distributed  Clutter 

Executive  Summary 

This  work  supports  the  Microwave  Radar  Branch’s  research  efforts  in  the  area  of  high 
grazing  angle  detection  of  small  maritime  targets  from  an  airborne  surveillance  platform.  It 
contributes  to  the  research  requirements  of  Task  07/040  (Support  to  AIR  70001).  Detection 
of  targets  is  a  major  function  of  radar  systems,  and  hence  detector  performance  is  an 
important  component  in  the  efforts  of  Task  07/040.  Detectors  are  designed  to  find  targets 
within  specific  clutter  models.  Recently  it  was  found  the  simple  Pareto  Distribution  fitted 
high  grazing  angle  clutter  quite  well.  This  was  validated  using  the  DSTO  Ingara  data, 
collected  in  a  trial  during  2004.  Hence  it  is  necessary  to  investigate  the  design  of  radar 
detection  schemes  in  such  a  clutter  environment. 

The  theory  of  spherically  invariant  random  processes  (SIRPs)  puts  this  problem  into  a 
formulation  that  allows  the  determination  of  detection  decision  rules.  The  Pareto  distri¬ 
bution  is  put  into  this  framework,  and  the  Neyman-Pearson  detector  is  specified.  Due  to 
intractibility  of  the  latter,  several  suboptimal  approximations  are  considered.  The  major 
contribution  of  this  work  is  to  examine  coherent  multilook  detection  under  the  assump¬ 
tion  of  a  Pareto  clutter  model.  Although  the  Pareto  SIRP  and  form  of  the  generalised 
likelihood  ratio  test  for  detection  has  been  presented  already  in  the  literature,  it  has  not 
been  applied  directly  to  the  multilook  problem  under  consideration. 

In  addition  to  this,  it  will  be  shown  that  for  Pareto  models  whose  parameters  are 
typical  of  clutter  obtained  from  a  radar  operating  in  vertical  polarisation,  the  whitening 
matched  filter  is  an  appropriate  suboptimal  detector.  This  is  useful  because  the  generalised 
likelihood  ratio  test  solution  depends  on  clutter  parameters,  unlike  the  whitening  matched 
filter. 


1  Future  Maritime  Patrol  and  Response  Capability. 
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Figures 

4.1  ROC  curves,  showing  the  effect  of  varying  the  Pareto  parameters.  Top  left 
plot  is  for  a  =  3  and  (3  =  1,  top  right  is  for  a  =  10  and  (3  =  1.  Bottom  left  is 

for  a  =  3  and  (3  =  10,  bottom  right  is  for  a  =  3  and  j3  =  100 .  12 

4.2  ROC  curves,  showing  the  WMF  becoming  a  better  suboptimal  decision  rule 
as  a  is  increased.  In  all  cases,  (3  =  0.01.  Top  left  is  for  a  =  3,  top  right  is  for 

a  =  10,  bottom  left  is  for  a  =  30  and  bottom  right  is  for  a  =  50 .  14 

4.3  ROC  curves  based  upon  the  Ingara  data  set  run34683.  Top  subplots  are  for 
vertical  polarisation,  bottom  two  are  for  horizontal  polarisation.  The  WMF 

is  suitable  only  for  the  vertically  polarised  case .  15 

4.4  ROC  curves  based  upon  Ingara  data  set  run34683,  with  azimuth  angle  190°. 

Top  two  subpolots  are  for  the  vertically  polarised  case,  while  bottom  two 
are  for  horizional  polarisation.  For  vertical  polarisation,  the  WMF  is  a  good 
approximation  to  the  GLRT.  For  the  horizontally  polarised  case,  increasing 

the  number  of  looks  seems  to  make  the  WMF  improve  in  performance.  ...  17 
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1  Preliminaries 


1.1  Introduction 

The  Pareto  distribution  [1,  2]  is  a  simple  power  law  model  that  has  been  proposed  as  an 
alternative  to  the  more  complex  high  resolution  radar  maritime  clutter  models  such  as  the 
K-  and  KK-Distributions  [3].  It  has  been  validated  both  for  the  case  of  low,  as  well  as 
high,  grazing  angle  sea  clutter  returns  [4,  5].  Its  advantage  is  that  it  is  a  much  simpler 
model  that  only  requires  two  parameters  to  be  fitted,  and  its  performance  matches  that 
of  other  modern  clutter  models.  Given  this,  it  is  thus  important  to  investigate  detection 
schemes  for  targets  in  Pareto  distributed  clutter.  In  order  to  achieve  this  it  is  necessary  to 
construct  the  multidimensional  clutter  and  signal  probability  distribution  functions  (pdfs). 
Such  an  exercise  is  quite  challenging  and  so  a  general  framework  is  required.  The  theory  of 
spherically  invariant  random  processes  (SIRPs)  [6-9]  provides  a  solution  to  this  problem. 
Such  processes  model  the  clutter  in  the  complex  domain  as  a  product  of  a  nonnegative 
random  variable  S,  and  a  zero  mean  complex  Gaussian  process  Q.  The  latter  process  and 
S  are  assumed  to  be  independent,  and  Q  has  a  N  x  N  covariance  matrix  £.  We  will  assume 
the  latter  is  semi-positive  definite,  so  that  its  inverse  possesses  a  Cholesky  decomposition, 
as  in  [10]. 

The  application  of  the  Pareto  distribution  to  sonar  detection  performance  is  outlined 
in  [7],  which  is  a  useful  reference  on  the  analysis  to  follow.  Here  we  focus  on  the  radar 
detection  problem,  based  upon  a  series  of  N  returns,  which  differs  significantly  from  the 
application  in  [7],  although  the  mathematical  derivations  of  detectors  is  similar. 

Coherent  multilook  radar  detection  for  targets  embedded  within  Pareto  distributed 
clutter  has  been  introduced  in  [11],  which  includes  generalised  likelihood  ratio  test  detec¬ 
tors  and  comparison  to  whitening  matched  filters.  The  Pareto  distribution  is  embedded 
within  a  spherically  invariant  random  process  (SIRP),  which  enables  the  decision  rules  to 
be  determined.  The  Pareto  distribution  has  been  shown  in  [12]  to  arise  as  the  intensity 
distribution  of  a  compound  Gaussian  model  for  which  a  linear  threshold  detector  is  opti¬ 
mal.  Other  recent  applications  of  the  Pareto  distribution  to  radar  include  [13,  14]  who  use 
the  fact  that  from  extreme  value  theory,  the  tail  of  any  distribution  can  be  modelled  as 
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a  generalised  Pareto  distribution.  Consequently,  CFAR  processes  are  derived  for  targets 
in  Pareto  distributed  clutter.  However,  the  resultant  threshold/false  alarm  probability 
relationship  is  dependent  on  both  Pareto  parameters,  which  is  a  significant  shortcom¬ 
ing.  A  Bayesian  approach  to  constructing  Pareto  CFAR  detectors  is  outlined  in  [15],  but 
the  resultant  threshold/false  alarm  probability  relationship  is  not  amenable  to  numerical 
methods.  Coherent  radar  detection  for  the  Pareto  distribution  has  also  been  considered 
in  [16],  from  the  point  of  view  of  a  compound  Gaussian  model  with  inverse  gamma  tex¬ 
ture1.  It  is  shown  that  the  CFAR  property  does  not  hold  in  general  for  coherent  detection 
schemes  considered.  However,  this  is  addressed  by  introducing  some  conditions  under 
which  a  CFAR  process  can  be  produced. 


1.2  Structure  and  Contributions 


The  report  does  not  present  a  detailed  analysis  of  SIRPs  but  instead  derives  the  ap¬ 
propriate  properties  directly  as  required.  This  is  also  specialised  to  the  problem  under 
consideration,  which  simplifies  the  form  of  the  detector  by  applying  a  whitening  approach 
first.  Section  2  formulates  the  multilook  detection  problem.  In  addition,  the  Pareto  SIRP 
is  introduced.  Although  this  is  presented  in  [7],  the  main  difference  is  the  performance  of 
multilook  detection  schemes  for  radar  are  considered  here.  The  optimal  Neyman-Pearson 
detector  is  constructed  in  Section  3  for  the  case  of  a  completely  known  target  model.  This 
assumption  is  then  relaxed,  by  constructing  the  generalised  likelihood  ratio  test  (GLRT),  as 
a  suboptimal  approximation  to  the  optimal  detector.  Section  4  analyses  the  performance 
of  the  GLRT  detector  by  comparing  it  with  a  whitening  matched  filter  approximation. 
Detector  performance  is  based  upon  a  Gaussian  target  model,  embedded  within  clutter 
whose  parameters  have  been  estimated  from  Ingara  data  sets.  This  is  in  contrast  to  the 
numerical  analysis  in  [16],  whose  Pareto  parameter  values  used  do  not  correspond  to  those 
estimated  from  real  clutter.  As  a  consequence  of  this,  it  will  be  shown  the  whitening 
matched  filter  can  be  used  as  a  good  suboptimal  detector  for  a  number  of  scenarios. 

xThis  paper  was  published  at  the  time  this  report  was  first  written 
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2  Detection  Problem  Structure 

2.1  Spherically  Invariant  Random  Processes 

Spherically  Invariant  Random  Processes  (SIRPs)  provide  a  general  formulation  of  the 
joint  density  of  a  non-Gaussian  random  process,  enabling  the  construction  of  densities  for 
Neyman-Pearson  detectors  [2,  17].  Traditionally,  a  SIRP  is  introduced  as  a  process  whose 
finite  order  subprocesses,  called  Spherically  Invariant  Random  Vectors  (SIRVs),  possess  a 
specific  density  [18,  19].  However,  due  to  an  equivalent  formulation,  we  can  specify  SIRVs 
and  SIRPs  in  a  manner  more  intuitive  to  the  modelling  of  radar  returns  as  follows.  Let 
c  =  {ci,  C2,  •  •  • ,  c/v}  be  the  complex  envelope  of  the  clutter  returns.  Then  this  vector  is 
called  SIRV  if  it  can  be  written  in  the  compound-Gaussian  formulation 

c  =  SG,  (2.1) 

where  the  process  Q  =  (Gu,  G2, . . . ,  Gn)T  is  a  zero  mean  complex  Gaussian  random  vec¬ 
tor,  or  multidimensional  complex  Gaussian  process,  and  S'  is  a  nonnegative  real  valued 
univariate  random  variable  with  density  fs  ■  The  latter  random  variable  is  assumed  to  be 
independent  of  the  former  process.  If  such  a  decomposition  exists  for  every  IV  £  IN  the 
complex  stochastic  process  {ci,C2,...}  is  called  spherically  invariant.  For  a  comprehen¬ 
sive  description  of  the  modelling  of  clutter  via  SIRPs  consult  [18].  What  is  clear  from  the 
formulation  (2.1)  is  that,  by  conditioning  on  the  random  variable  S,  we  can  write  down 
the  density  of  c  as  a  convolution.  This  will  be  shown  explicitly  in  the  analysis  to  follow. 

2.2  General  Formulation 

The  detection  decision  problem  is  now  formulated.  We  assume  that  we  have  a  SIRP  model 
for  the  clutter  c  as  specified  through  the  product  formulation  (2.1).  Suppose  the  radar 
return  is  z,  which  is  a  complex  N  x  1  vector.  Then  the  coherent  multilook  detection 
problem  can  be  cast  in  the  form 

Hq  :  z  =  c  against  H\  :  z  =  Rp  +  c,  (2.2) 

where  all  complex  vectors  are  N  x  1,  and  Hq  is  the  null  hypothesis  (return  is  just  clutter) 
and  H\  is  the  alternative  hypothesis  (return  is  a  mixture  of  signal  and  clutter).  Statistical 
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hypothesis  testing  is  outlined  in  [2].  Here,  the  vector  p  is  the  Doppler  steering  vector, 
whose  components  are  given  by  p(j)  =  e--727™-'0,  for  j  £  {0, 1, . . .  ,  N  —  1},  where  fr>  is  the 
normalised  target  Doppler  frequency,  and  —0.5  <  fr>  <  0.5.  It  will  be  assumed  that  this 
is  completely  known.  The  complex  random  variable  R  accounts  for  target  characteristics, 
and  \R\  is  the  target  amplitude.  Suppose  the  zero  mean  Gaussian  process  has  a  covariance 
matrix  cov(G)  =  1E(GGH)  =  S.  Since  it  is  a  covariance  matrix,  its  inverse  will  exist.  It 
can  be  shown  that  the  inverse  of  £  is  symmetric  and  positive  definite.  If  we  assume  it 
is  semi-definite  positive,  then  we  can  apply  a  whitening  filter  approach  to  simplify  the 
detection  problem.  Hence  we  suppose  the  Cholesky  Factorisation  exists  for  £_1,  so  that 
there  exists  a  matrix  A  such  that  £_1  =  AH A. 

A  SIRP  is  unaffected  by  a  linear  transformation  [6].  This  means  that  applying  a 
linear  operator  to  the  clutter  process  alters  the  complex  Gaussian  component,  but  the 
characteristic  function  of  the  SIRP  is  preserved.  Hence  in  the  literature,  a  whitening 
approach  is  often  applied  to  the  detection  problem  of  interest.  Note  that,  by  applying  the 
Cholesky  factor  matrix  A  to  the  statistical  test,  we  can  reformulate  (2.2)  in  the  statistically 
equivalent  form 

Ho  :  r  =  n  against  H\  :  r  =  Ru  +  n,  (2.3) 

where  r  =  Az ,  n  =  Ac  and  u  =  Ap. 

The  transformed  clutter  process  n  =  SAG,  and  AG  is  still  a  multidimensional  complex 
Gaussian  process,  with  zero  mean  but  covariance  cov(AG)  =  E (AQGH AH)  =  AY*AH . 

Since  £_1  =  AH A,  it  follows  that  Inxn  =  ^AH A,  from  which  it  is  not  difficult  to 
deduce  that  A  =  ( ATiAh)A .  Let  B  =  ATiAh ,  then  note  that  B2  =  B,  so  that  B  is  idem- 
potent.  Also,  it  follows  that  B  must  be  invertible,  since  det(B)  =  det(A)det(T,)det(AH ) 
and  det{Ti)  /  0  and  det(£-1)  =  det(AH)det(A)  /  0.  Thus  it  follows  that  B  must  be 
the  N  x  N  identity  matrix.  Consequently,  we  conclude  that  the  transformed  Gaussian 
process  AG  has  zero  mean  and  covariance  matrix  this  identity  matrix.  Hence  the  clutter, 
conditioned  on  the  variable  S,  is  completely  decorrelated  through  this  linear  transform. 

We  write  this  as  AG  =  CN( 0,  Inxn)-  Observe  that  n\S  is  still  complex  Gaussian  with 
zero  mean  vector  but  covariance  matrix  s2InxNi  hence  its  density  is 

fn\S=s(x )  =  nNg2Ne  ^ 
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where  the  norm  in  (2.4)  is  the  complex  Euclidean  norm  on  vectors.  Hence  by  using 
conditional  probability  [20] 


fn(x)  = 


fn\S=s(x)fs(s)ds 


I  o  1TNS2N 


1  e-s'2^2  fs(s)ds. 


Hence  if  we  define  a  function  h^ip)  by 

/•OO 

hN{p)  =  /  s~2Ne~s  2pfs(s)ds, 

Jo 

then  the  clutter  joint  density  can  be  written  in  the  compact  form 

fnix)  =  -^hN(\\x\\2). 


(2.5) 


(2.6) 


(2.7) 


The  function  h jv  defined  in  (2.6)  is  of  paramount  importance  in  SIRP  theory,  as  all  densities 
of  interest  are  expressed  in  terms  of  it.  It  is  thus  called  the  characteristic  function.  Much 
of  the  literature  is  devoted  to  determining  Hn  and  fs  pairs  for  particular  desired  clutter 
models  [6,  19]. 


Next  we  derive  the  marginal  amplitude  and  intensity  distributions  of  the  complex 
clutter  process  n,  which  are  intimately  related  to  the  special  function  /ijv-  Let  the  fcth 
element  of  n  be  n&,  so  that  =  SAG where  Gk  is  the  fcth  component  of  G ,  for  1  <  k  <  N. 
Then  the  complex  Gaussian  density  of  AGk  is 

IasM  =  (2-8) 

Hence,  equivalently,  AGk  is  a  complex  Gaussian  process  with  mean  zero  and  covariance 
matrix  -7/2x2-  Consequently,  it  follows  that  its  amplitude  \AQk\  has  a  Rayleigh  distribution 
with  parameter  A=.  Hence,  using  the  fact  that  S  and  AGk  are  independent,  and  the 
definition  of  the  Rayleigh  density,  we  can  show  that 

f°°  2 1  t2 

/|nfc|(t)=/  —e~^fs{s)d.s  =  2thi{t2),  (2.9) 

Jo  s 

yielding  the  marginal  amplitude  distribution  for  each  k.  Transforming  to  the  intensity 
domain,  it  is  not  difficult  to  show 

f\nk\2(t)  =  f\nk\{Vt)  =  hi(t ),  (2.10) 
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also  for  each  k. 

The  importance  of  the  results  (2.9)  and  (2.10)  is  that  they  indicate  how  an  arbitrary 
distribution  can  be  embedded  within  a  SIRP  model.  The  key  to  the  determination  of  a 
relevant  SIRP  for  a  given  detection  problem  is  the  specification  of  the  random  variable 
S  through  its  density  fs •  For  a  desired  marginal  distribution,  we  can  form  an  integral 
equation  using  (2.10)  and  (2.6)  with  N  =  1.  The  literature  contains  extensive  discussion 
of  this  problem;  in  particular,  [6]  contains  guidelines  on  performing  this,  as  well  as  case 
studies  for  most  desired  marginal  distributions. 


2.3  Pareto  SIRP 


The  Pareto  SIRP  is  defined  to  be  that  whose  marginal  intensity  distributions  are  Pareto, 


with  given  parameters  a  and  /3.  Recall  that  X  =  Pa(a,  (3)  if  its  density  is  given  by 

0  t  <  p, 


fx(t )  = 


t>0. 


(2.11) 


The  parameter  a  is  called  the  distribution’s  shape  parameter,  while  /3  is  called  the  scale 
parameter.  The  latter  specifies  where  on  the  real  line  the  distribution’s  support  begins, 
which  is  the  interval  [/3,  oo).  The  shape  parameter  controls  how  fast  the  distribution’s  tail 
decreases.  In  order  to  specify  a  multidimensional  process  of  the  form  (2.1),  we  need  to 
specify  an  appropriate  density  of  S.  The  choice  below  has  been  based  upon  the  results  in 

[7]- 

Consider  the  function  fs(s)  defined  by 


2/3"  _2a—l  - 0s~ 2 
fs{s)  =  -p--s  2q  le  Ps  , 

F(«) 


(2.12) 


for  s  >  0,  where  a  and  /3  are  positive.  An  application  of  L’Hopital’s  Rule  shows 
lims_>0+  fs(s)  =  0,  so  we  can  set  /$■(())  equal  to  zero  to  make  fs  continuous  on  the  half 
interval  [0,  oo).  Clearly  fs{s)  >  0  Vs  >  0  and  by  a  change  of  variables  to  u  =  /3s-2, 

r  oo  /3a~  i 


/  fs(s)ds  =  - 

/o  r(«)  Jo 


°o  /u\  a-l 


P 


e  Udu  =  1, 


(2.13) 


using  the  definition  of  the  Gamma  Function.  Hence  the  function  fs{s)  is  a  valid  density 
on  [0,  oo). 
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It  is  not  difficult  to  show  that  its  first  and  second  moments  are  given  by 


E (S)  =  ^T(f;  -  °'5)  and  E(S2)  =  '3 


F(a) 


a  —  1 


(2.14) 


Using  this  density,  we  determine  the  characteristic  function  hiy .  By  applying  the 
transform  u  =  [p  +  /3]s-2,  we  can  show 


hN{jp)  = 


2/3 ° 

I>) 


-2N-2a-le-\p+/3]s  2 


oa  /*oo 

=  _ - _  /  uN+a~le~udu 

T(a)\p  +  (3]^J0 

(3°T(N  +  a) 

r(a)[p  +  /3]w+«’ 

where  the  definition  of  the  Gamma  Function  has  been  applied  as  previously. 


(2.15) 


Next  we  construct  the  marginal  intensity  distribution.  By  an  application  of  (2.10),  we 
see  that  by  setting  IV  =  1  in  (2.15), 

/KI’W  =  | r^TT>  (2.16) 

showing  the  marginal  distributions  corresponding  to  the  whitened  SIRP  with  the  choice 
of  fs(s )  in  (2.12)  being  shifted  Pareto  distributions.  Note  that  the  linear  transform  z k  = 
|^-fc|2  +  (3  yields  a  density 

-  m  =  <2-17) 

for  t  >  f3,  which  is  the  standard  form  of  the  Pareto  density.  Hence,  the  clutter  process 
is  modelled  as  a  product  of  a  random  variable  generated  by  S  with  density  (2.12),  and  a 
complex  Gaussian  process  with  covariance  matrix  £.  Note  that  since  the  Gaussian  process 
is  zero  mean,  cov(e)  =  E(S’2)S,  implying  we  can  produce  a  desired  covariance  structure 
for  our  clutter  process  by  scaling  the  Gaussian  process  by  the  mean  squared  value  of  S. 


3  Neyman-Pearson  Detectors 

We  now  turn  to  the  determination  of  optimal  and  suboptimal  decision  rules  using  the 
Neyman-Pearson  Lemma  [17]  and  the  SIRP  structure  derived  in  the  previous  section.  We 
begin  with  the  case  of  a  fully  specified  target  model. 
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3.1  Case  of  a  Completely  Known  Target 


It  is  useful  to  begin  with  the  case  of  a  fixed  known  target  model,  which  means  we  assume 
that  R  is  a  fixed  constant  that  is  completely  specified. 


Under  Hi,  r  =  Ru  +  n,  and  conditioned  on  S,  this  is  still  complex  Gaussian,  with  the 
same  covariance  as  n  (since  Ru  is  constant),  but  its  mean  is  shifted  by  the  vector  Ru. 
Hence  the  density  under  H\  is 

fair)  =  -^hN{\\r  -  Ru\\2).  (3.1) 

7 TJV 


Consequently,  the  Neyman-Pearson  optimal  detector  is  the  ratio  of  the  densities  under  H\ 
and  Hq.  The  density  under  H$  is  given  by  (2.7),  with  an  application  of  (2.15).  Hence  the 
likelihood  is 


L(r)  = 


+  P\ 


N-\-a 


J|r  —  Ru\\2  +  (3]N+a '' 
using  (3.1).  Thus,  the  likelihood  ratio  test  takes  the  form 


M(r)  = 


+  (3 


S1 

t.  T 


(3.2) 


(3.3) 


|r  —  .Rtt||2  +  (3  h0 

where  r  is  the  detection  threshold,  which  can  be  determined  numerically  from  the  false 


alarm  specification.  The  notation  X  ^  Y  means  that  we  reject  H$  if  and  only  if  X  >  Y. 

H0 

In  all  reasonable  applications,  we  will  not  have  knowledge  of  R,  and  so  we  must  extend 
this  result  to  account  for  an  unknown  R.  This  is  the  subject  of  the  next  subsection. 


3.2  Unknown  Target:  Generalised  Likelihood  Ratio  Test 

In  this  case  we  assume  R  is  unknown  but  constant  within  a  given  range  profile.  This  means 
we  can  apply  the  GLRT  to  first  estimate  this  parameter,  and  then  apply  the  estimate  to 
the  test  (3.3).  The  methodology  of  GLRT  is  described  in  [7,  10].  We  essentially  produce 
a  suboptimal  detector  for  the  case  of  an  unknown  target,  based  upon  the  optimal  decision 
rule  formulated  for  a  completely  known  target  model. 

We  call  this  estimate  i?,  which  is  taken  to  be  the  maximum  likelihood  estimator  [2] .  It 
is  chosen  to  minimise  the  quantity  ||r  —  Ru\\2.  This  is  because,  in  view  of  the  density  (3.1) 
and  the  definition  of  h]\r  in  (2.6),  maximising  the  density  with  respect  to  R  is  equivalent 
to  minimisation  of  this  quantity. 
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By  applying  a  simple  expansion,  we  can  show 

|| r  -  Ru\\2  =  ||r||2  +  |i?|2||tt||2  -  2 $t(Runr),  (3.4) 


where  3?(z)  is  the  real  part  of  the  complex  number  z.  Recalling  the  proof  of  the  Cauchy- 
Schwarz  inequality  (see  [21]  for  example),  this  is  minimised  with  the  choice  R  =  jjUjp. 
Substituting  this  in  (3.4)  shows  that 


„  \vHr\ 2 

min  ||r  —  Ru\\2  =  Mr*  —  Ru||2  =  ||r||2 - - — 4-. 


R 


An  application  of  (3.5)  to  (3.3)  results  in  the  GLRT 


L(r)  = 


+  (3 


§ 

I  no  \uHr\2  .  n 

Ml  +  0 


r. 


(3.5) 


(3.6) 


This  result  is  also  derived  in  [7]  in  the  context  of  sonar  signal  processing.  However,  our 
application  is  significantly  different  from  that  in  the  latter.  The  decision  rule  (3.6)  gives 
a  suboptimal  detector  for  the  detection  of  a  target  in  clutter  that  has  Pareto  marginal 
distributions.  The  next  Section  will  analyse  its  perfomance. 


4  Detector  Performance  Analysis 

The  performance  of  the  GLRT  detector  (3.6)  is  now  examined.  It  will  also  be  compared  to 
the  performance  of  the  whitening  matched  filter  (WMF),  which  can  serve  as  a  suboptimal 
decision  rule.  Determining  when  the  WMF  performs  well  is  important  because  it  does 
not  depend  on  the  clutter  parameters,  as  does  the  GLRT  solution.  Four  examples  will 
be  considered:  the  first  two  explore  the  role  of  the  Pareto  clutter  parameters  on  detector 
performance.  The  second  two  examine  the  performance  of  the  decision  rules  when  applied 
to  scenarios  with  Pareto  parameters  estimated  from  the  2004  Ingara  data  sets  [22],  How¬ 
ever,  the  clutter  is  simulated  using  the  estimated  parameters,  to  determine  the  receiver 
operating  characteristic  curves.  Details  of  the  Ingara  clutter  are  not  repeated  here,  but 
are  available  in  a  number  of  references.  Refer  to  [3]  and  [22]  for  specific  details  of  the  data 
and  trial,  and  [5]  for  information  on  the  Pareto  fit  to  the  Ingara  data. 

For  a  given  false  alarm  probability,  the  detection  threshold  r  has  been  estimated  using 
Monte  Carlo  simulation.  The  number  of  independent  and  identically  distributed  samples 
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used  is  106  in  all  cases.  For  a  given  SCR,  the  detection  probability  has  also  been  estimated 
using  Monte  Carlo  simulation,  where  for  each  given  SCR,  106  independent  and  identically 
distributed  samples  have  been  used  to  estimate  the  detection  probability. 

4.1  Signal  to  Clutter  Ratio 

The  signal  to  clutter  ratio  (SCR)  measures  the  signal  power  relative  to  the  clutter  power. 
It  is  defined  as  SCR  =  S1/S2,  where 

5i  =  E|R«|2  =  E[(i?«)H(R«)]  =  E|R|2||«||2.  (4.1) 

Since  R  is  assumed  to  be  bivariate  Gaussian,  it  follows  that  \R\  has  a  Rayleigh  distribution 
with  parameter  a,  and  so 

E|R|2  =  var(|R|)  +  (E(|i?|))2  =  2u2.  (4.2) 

Next,  the  clutter  has  mean  squared  value 

S2  =  E|n|2  =  E(,S2)E|A£|2.  (4.3) 

But  recall  the  effect  of  the  matrix  A  is  that  it  decorrelates  the  Gaussian  component  Q . 
Hence  it  is  simple  to  show  E|A£|2  =  N,  and  so  S2  =  A^E(S'2).  Consequently  it  follows 
that 

_  2fr2||u||2  _  2a2pHT,~1p  _  2<j2pHY,~1p(a  -  1) 

E  (S2)N  E  (S2)N  (3  N  ’  [  ) 

using  (2.14).  In  order  to  assess  detector  performance,  the  SCR  is  varied  from  -10  to  +30 
dB,  and  a  Gaussian  target  model  is  used,  whose  parameters  are  determined  from  (4.4). 

4.2  Whitening  Matched  Filter 

As  a  comparison,  the  whitening  matched  filter,  or  Gaussian  optimal  detector  is  used  as  a 
suboptimal  approximation.  For  a  return  r,  this  decision  rule  takes  the  form 

\uHr\2  ^  rwmf,  (4.5) 

H0 

where  Twmj  is  the  detector’s  threshold. 
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4.3  Receiver  Operating  Characteristics 

A  series  of  receiver  operating  characteristic  (ROC)  curves  are  now  provided  to  assess 
detector  performance.  Throughout,  a  Toeplitz  structure  will  be  assumed  for  the  clutter 
covariance  matrix  £.  In  particular,  it  will  be  assumed  that  T,(j,k)  =  for  each 

j,  k  €  {0, 1,  2, . . .  N  —  1},  and  0  <  k  <  1.  As  pointed  out  in  [3],  the  azimuth  angle  of  Ingara 
is  the  angle  the  radar  looks  towards.  For  the  Ingara  trial,  the  wind  speed  is  reported  to 
have  been  7.1  metres  per  second,  in  a  direction  of  47°,  for  the  Australian  Government’s 
Bureau  of  Meterology  estimates.  This  means  that  the  upwind  angle  is  227°,  which  is  the 
point  where  the  clutter  will  be  strongest  [3]. 

All  ROC  curves  show  the  performance  of  the  suboptimal  detectors  relative  to  the 
optimal  decision  rule  (3.3),  denoted  OPT  throughout. 


4.3.1  Example  1 

The  first  set  of  ROC  curves  we  examine  explore  the  effect  of  the  Pareto  parameters  on 
detector  performance.  Figure  4.1  is  for  the  case  where  the  number  of  looks  is  N  =  5,  the 
normalised  Doppler  frequency  is  fo  =  0.5,  the  probability  of  false  alarm  is  10-6  and  the 
Toeplitz  factor  is  n  =  0.1.  The  top  two  subplots  in  figure  4.1  show  the  effect  of  increasing 
the  Pareto  shape  parameter  a.  The  top  left  subplot  is  for  the  case  where  a  =  3  and  (3=1. 
The  top  right  subplot  is  for  the  case  where  a  =  10  and  (3=1.  These  two  plots  show 
a  phenomenon  observed  quite  frequently  in  the  ROC  analysis.  As  the  shape  parameter 
is  increased,  the  WMF  becomes  a  good  approximation  to  the  GLRT  detector.  This  is 
probably  because  the  clutter  is  becoming  more  Gaussian  like.  The  bottom  two  subplots 
in  figure  4.1  show  the  effect  of  changing  the  scale  parameter,  while  the  shape  parameter 
is  held  constant.  The  bottom  left  subplot  is  for  a  =  3  and  (3  =  10,  while  the  bottom  right 
subplot  is  for  a  =  3  and  (3  =  100.  This  was  observed  in  other  simulations:  it  appears  that 
the  scale  parameter  does  not  affect  the  detector’s  performance  significantly.  We  observe 
that  the  GLRT  has  good  performance  throughout. 
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ROC  Curve:  Pareto  Multilook  Coherent  Detection,  N  =  5  and  fD  =  0.5 


SCR  (dB) 


ROC  Curve:  Pareto  Multilook  Coherent  Detection,  N  =  5  and  fQ  =  0.5 


ROC  Curve:  Pareto  Multilook  Coherent  Detection,  N  =  5  and  fD  =  0.5 


SCR  (dB) 


ROC  Curve:  Pareto  Multilook  Coherent  Detection,  N  =  5  and  fD  =  0.5 


Figure  f.l:  ROC  curves,  showing  the  effect  of  varying  the  Pareto  parameters.  Top  left 
plot  is  for  a  =  3  and  (3  =  1,  top  right  is  for  a  =  10  and  (3  =  1.  Bottom  left  is  for  a  =  3 
and  (3  =  10,  bottom  right  is  for  a  =  3  and  (3  =  100. 
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4.3.2  Example  2 

As  a  second  example,  the  effect  of  varying  a,  while  the  scale  parameter  is  held  fixed, 
is  examined.  Figure  4.2  is  for  the  case  where  N  =  10,  fo  =  1,  false  alarm  probability 
1CD6  and  Toeplitz  factor  k  =  0.9  is  considered.  The  scale  parameter  is  fixed  at  (3  =  0.01, 
while  the  shape  parameter  is  increased.  This  type  of  scale  parameter  has  been  found  to 
be  typical  for  horiziontally  polarised  Ingara  clutter  returns,  when  modelled  by  a  Pareto 
distribution.  The  top  left  plot  in  figure  4.2  is  for  a  =  3,  the  top  right  is  for  a  =  10,  the 
bottom  left  is  for  a  =  30  and  the  bottom  right  subplot  is  for  a  =  50.  This  example  shows 
the  WMF  becoming  a  very  good  suboptimal  decision  rule  as  a  increases.  For  large  a,  it 
matches  the  GLRT  very  closely.  As  before,  the  GLRT  is  always  a  good  approximation  to 
the  optimal  decision  rule. 


4.3.3  Example  3 

The  third  example  considered  uses  Pareto  parameters  estimated  from  the  Ingara  data  sets. 
In  this  case,  the  parameters  are  sourced  from  run34683,  at  an  azimuth  angle  of  225°.  This 
is  only  2°  off  from  the  upwind  direction,  and  so  the  clutter  should  be  strong.  Throughout, 
the  normalised  Doppler  frequency  has  been  set  to  fo  =  0.4,  the  Toeplitz  factor  k  =  0.5 
and  false  alarm  probability  of  10~6. 

The  top  two  subplots  of  figure  4.3  are  for  the  vertically  polarised  case,  where  the  Pareto 
parameters  have  been  estimated  to  be  a  =  11.3930  and  (3  =  0.3440.  The  top  left  subplot 
has  number  of  looks  N  =  2,  while  the  top  right  is  for  N  =  8.  In  this  case,  the  WMF  is  a 
good  approximation  to  the  GLRT. 

The  bottom  two  subplots  are  for  the  horizontally  polarised  case,  where  the  estimated 
Pareto  parameters  are  a  =  4.7241  and  (3  =  0.0466.  The  bottom  left  subplot  is  for  N  =  2 
and  the  bottom  right  subplot  is  for  N  =  8.  Here  the  WMF  is  not  such  a  good  approx¬ 
imation.  In  many  of  the  Ingara  data  sets  examined,  this  was  found  to  be  the  case:  the 
vertically  polarised  case  resulted  in  the  WMF  being  a  good  suboptimal  detecor,  while  for 
horizontally  polarised  clutter,  the  GLRT  is  more  suitable.  The  GLRT  tends  to  only  incur 
a  small  detection  loss  throughout. 
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ROC  Curve:  Pareto  Multilook  Coherent  Detection,  N  =  10  and  fD  =  1 


SCR  (dB) 


ROC  Curve:  Pareto  Multilook  Coherent  Detection,  N  =  10  and  fQ  =  1 


ROC  Curve:  Pareto  Multilook  Coherent  Detection,  N  =  10  and  fD  =  1 


SCR  (dB) 


ROC  Curve:  Pareto  Multilook  Coherent  Detection,  N  =  10  and  fD  =  1 


Figure  f..2:  ROC  curves,  showing  the  WMF  becoming  a  better  suboptimal  decision  rule 
as  a  is  increased.  In  all  cases,  (3  =  0.01.  Top  left  is  for  a  =  3,  top  right  is  for  a  =  10, 
bottom  left  is  for  a  =  30  and  bottom  right  is  for  a  =  50. 
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ROC  Curve:  Pareto  Multilook  Coherent  Detection,  N  =  2  and  fD  =  0.5 


SCR  (dB) 


ROC  Curve:  Pareto  Multilook  Coherent  Detection,  N  =  8  and  fQ  =  0.5 
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SCR  (dB) 


ROC  Curve:  Pareto  Multilook  Coherent  Detection,  N  =  8  and  fD  =  0.5 


Figure  4-3:  ROC  curves  based  upon  the  Ingara  data  set  run34683.  Top  subplots  are  for 
vertical  polarisation,  bottom  two  are  for  horizontal  polarisation.  The  WMF  is  suitable  only 


for  the  vertically  polarised  case. 
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4.3.4  Example  4 

As  a  final  example,  the  ROC  performance  of  detectors  with  Pareto  parameters  sourced 
from  run34683,  at  an  azimuth  of  190°,  is  considered.  In  this  case,  the  clutter  should 
be  weaker  because  the  radar  is  focused  in  the  downwind  direction  (37°  discrepancy  from 
downwind).  In  this  example,  the  normalised  Doppler  frequency  has  been  set  to  fo  =  0.5, 
the  Toeplitz  factor  is  k  =  0.01  and  false  alarm  probability  of  1CU6.  For  the  vertically 
polarised  case,  a  =  8.9028  and  (3  =  0.1109.  For  the  horizontally  polarised  case,  a  =  4.9031 
and  (3  =  0.0209.  Figure  4.4  shows  the  results.  The  top  two  plots  are  for  the  vertically 
polarised  case,  while  the  bottom  two  are  for  the  horizontal  polarisation.  The  number  of 
looks  is  varied  from  N  =  2  to  N  =  8,  moving  left  to  right  as  illustrated.  For  the  vertically 
polarised  case,  we  see  that  the  WMF  is  a  good  approximation,  while  for  the  horizontally 
polarised  case,  it  becomes  better  as  the  number  of  looks  increases.  Observe  the  GLRT  has 
consistent  performance. 


4.4  Detector  Computational  Load  Considerations 

Due  to  the  fact  that  a  radar  system  must  process  many  sets  of  returns  rapidly,  the  rel¬ 
ative  computation  loads  of  detectors  is  an  important  consideration  in  detector  design. 
The  WMF  detector  (equation  4.5)  requires  three  operations  (multiplication  of  Hermitian 
transpose  vector  and  whitened  return,  followed  by  modulus  and  squaring).  It  is  also  in¬ 
dependent  of  the  clutter  parameters.  By  contrast,  the  GLRT  has  a  higher  computational 
load  (see  equation  3.6).  Given  a  whitened  return  r,  the  GLRT  involves  calculation  of  the 
WMF,  normalisation  of  it  and  then  passing  the  result  into  a  ratio  function  which  requires 
computation  of  the  modulus  squared  of  the  return  r.  This  can  be  done  in  a  total  of  about 
6  operations.  However,  these  are  not  computationally  intense  operations? 

As  noted  previously,  the  performance  of  the  GLRT  was  matched  by  that  of  the  WMF 
in  a  number  of  cases  considered.  Given  this,  the  WMF  could  be  used  as  a  computationally 
cheaper  solution. 

2Examples  of  a  high  computational  load  is  evaluation  of  Bessel  functions  and  other  integrals,  such  as 
is  experienced  with  similar  detectors  for  the  case  of  K-  and  KK-Distributed  clutter. 
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ROC  Curve:  Pareto  Multilook  Coherent  Detection,  N  =  2  and  fD  =  0.5  ROC  Curve:  Pareto  Multilook  Coherent  Detection,  N  =  8  and  fD  =  0.5 


ROC  Curve:  Pareto  Multilook  Coherent  Detection,  N  =  2  and  fD  =  0.5 


SCR  (dB) 


ROC  Curve:  Pareto  Multilook  Coherent  Detection,  N  =  8  and  fD  =  0.5 


Figure  4-4:  ROC  curves  based  upon  Ingara  data  set  run34683,  with  azimuth  angle  190°. 
Top  two  subpolots  are  for  the  vertically  polarised  case,  while  bottom  two  are  for  horizional 
polarisation.  For  vertical  polarisation,  the  WMF  is  a  good  approximation  to  the  GLRT. 
For  the  horizontally  polarised  case,  increasing  the  number  of  looks  seems  to  make  the  WMF 
improve  in  performance. 
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5  Conclusions 

This  report  examined  coherent  multilook  detection  for  targets  in  Pareto  distributed  clutter. 
By  introducing  the  Pareto  SIRP,  the  Neyman-Pearson  Lemma  yielded  the  theoretical 
optimal  decision  rule.  From  this,  the  suboptimal  GLRT  was  constructed.  It  was  shown 
to  have  small  detection  loss  in  a  series  of  examples,  some  of  which  have  been  based  upon 
Pareto  clutter  estimates  sourced  from  Ingara  data  sets.  The  WMF  was  observed  to  be  a 
good  suboptimal  approximation  in  a  number  of  simulations,  suggesting  it  can  be  used  as  a 
computationally  simpler  solution.  The  relationship  between  the  WMF  and  GLRT  will  be 
explored  in  subsequent  work,  hopefully  providing  some  general  guidelines  when  the  WMF 
is  comparable  to  the  GLRT.  This  will  be  advantageous  from  the  perspective  that  it  does 
not  depend  on  the  Pareto  scale  parameter,  as  does  the  GLRT. 
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